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Gluon jets with a mean energy of 22 GeV and purity of 95% are selected from hadronic Z° 
decay events produced in e + e~ annihilations. A subsample of these jets is identified which 
exhibits a large gap in the rapidity distribution of particles within the jet. After imposing 
the requirement of a rapidity gap, the gluon jet purity is 86%. These jets are observed to 
demonstrate a high degree of sensitivity to the presence of color reconnection, i.e. higher order 
QCD processes affecting the underlying color structure. We use our data to test three QCD 
models which include a simulation of color reconnection: one in the Ariadne Monte Carlo, one 
in the Herwig Monte Carlo, and the other by Rathsman in the Pythia Monte Carlo. We find 
the Rathsman and Ariadne color reconnection models can describe our gluon jet measurements 
only if very large values are used for the cutoff parameters which serve to terminate the parton 
showers, and that the description of inclusive Z° data is significantly degraded in this case. 
We conclude that color reconnection as implemented by these two models is disfavored. The 
signal from the Herwig color reconnection model is less clear and we do not obtain a definite 
conclusion concerning this model. In a separate study, we follow recent theoretical suggestions 
and search for glueball-like objects in the leading part of the gluon jets. No clear evidence is 
observed for these objects. 
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1 Introduction 



Rapidity y, denned by y = \ In J with E the energy of a particle and p\\ the component 

of its 3- momentum along an axis 1 , is one of the most common variables used to characterize 
the phase space distribution of particles in high energy collisions. Of current interest (see for 
example [2]) are events with a so-called rapidity gap, namely events in which two populated 
regions in rapidity are separated by an interval devoid of particles. High energy collisions are 
often characterized by the formation of quark and gluon jets, i.e. collimated streams of hadrons 
associated with the hard scattering of quarks and gluons, respectively. Most recent interest in 
rapidity gaps has focused on a class of events in electron-proton [3] and proton-antiproton [4] 
collisions with large rapidity gaps between jets: these events are interpreted as arising from 
the exchange of a strongly interacting color singlet object, such as a pomeron [5], between the 
underlying partonic constituents of the event. 

Another source of rapidity gaps is color reconnection (CR), i.e. a rearrangement of the 
underlying color structure of an event from its simplest configuration, in which a color flux tube 
or "string" is stretched from a quark to an antiquark through intermediate gluons in a manner 
such that string segments do not cross (a so-called planar diagram, see Fig. la), to a more 
complex pattern in which some segments can either cross or else appear as disconnected entities 
whose endpoints are gluons (Fig. lb). Diagrams with color reconnection represent higher order 
processes in Quantum Chromodynamics (QCD), suppressed by order 1/Nq compared to planar 
diagrams, where Nq = 3 is the number of colors. In models of hadron production such as the 
Lund string model [6], the flux tubes hadronize. In events with a disconnected gluonic string 
segment as in Fig. lb, a rapidity gap can form between the isolated segment - often the leading 
(highest rapidity) part of a gluon jet - and the rest of the event. Thus rapidity gaps in gluon 
jets can provide a sensitive means to search for effects of color reconnection. Color reconnection 
has been a topic of considerable recent interest because of its potential effects in fully hadronic 
decays of W + W~ events produced in electron-positron (e + e~) collisions [7], introducing an 
uncertainty in the measurement of the W boson mass at LEP [8]. 

Recently [9], gluon jets with a rapidity gap were also proposed as a potentially favorable 
environment for the production of color singlet bound states of gluons, such as glueballs, through 
diagrams like Fig. lb in which the isolated gluonic system represents a hadronic resonance. 

Previous studies of rapidity gaps in e + e~ hadronic annihilations were based on inclusive Z° 
events and separated two- and three-jet events from Z° decays [10]. The rapidity distribution 
of charged particles in gluon jets was used to test models of color reconnection in [11]. There 
are no previously published experimental studies on gluon jets with a rapidity gap. 

In this paper, we study gluon jets with rapidity gaps, produced in three-jet quark- ant iquark- 
gluon (qqg) events from e + e~ hadronic Z° decays. The gluon jets are identified through "anti- 
tagging," using displaced secondary vertices from B hadrons to identify the quark and anti- 
quark jets. The data were collected with the OPAL detector at the LEP e + e~ storage ring at 
CERN. We measure the charged particle multiplicity, total electric charge, and distributions of 
invariant mass in the leading part of the gluon jets. 

Usually the thrust [1], jet, or beam axis. 
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Figure 1: Schematic illustrations of events with (a) standard "planar" color flow and (b) re- 
connection. The hatched regions represent color flux tubes or "strings" stretched between the 
quark q, antiquark q and gluons g. 

2 Detector and data sample 



The OPAL detector is described in detail elsewhere [12, 13]. OPAL operated from 1989 to 
2000 and subsequently was dismantled. The tracking system consisted of a silicon microvertex 
detector, an inner vertex chamber, a large volume jet chamber, and specialized chambers at the 
outer radius of the jet chamber to improve the measurements in the z-direction. 2 The tracking 
system covered the region | cos#| < 0.98 and was enclosed by a solenoidal magnet coil with an 
axial field of 0.435 T. Electromagnetic energy was measured by a lead-glass calorimeter located 
outside the magnet coil, which also covered | cos6>| <0.98. 

The present analysis is based on a sample of about 2 722 000 hadronic annihilation events, 
corresponding to the OPAL sample collected within 3 GeV of the Z° peak from 1993 to 1995. 
This sample includes readout of both the r-0 and z coordinates of the silicon strip microvertex 
detector [13]. The procedures for identifying hadronic annihilation events are described in [14]. 

We employ the tracks of charged particles reconstructed in the tracking chambers and 
clusters of energy deposited in the electromagnetic calorimeter. Charged tracks are required to 
have at least 20 measured points (of 159 possible) in the jet chamber, or at least 50% of the 
number of points expected based on the track's polar angle, whichever is larger. In addition, 
the tracks are required to have a momentum component perpendicular to the beam axis greater 
than 0.05 GeV/c, to lie in the region | cos#| < 0.96, to point to the origin to within 5 cm in the 
r-0 plane and 30 cm in the z direction, and to yield a reasonable x 2 P er degree-of-freedom for 

2 Our right handed coordinate system is defined so that z is parallel to the e~ beam axis, x points towards 
the center of the LEP ring, r is the coordinate normal to the beam axis, (f> is the azimuthal angle around the 
beam axis with respect to x, and 9 is the polar angle with respect to z. 
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the track fit in the r-0 plane. Electromagnetic clusters are required to have an energy greater 
than 0.10 GeV if they are in the barrel section of the detector (| cos#| <0.82) or 0.25 GeV if 
they are in the endcap section (0.82 < | cos#| < 0.98). A matching algorithm [15] is employed to 
reduce double counting of energy in cases where charged tracks point towards electromagnetic 
clusters. Specifically, if a charged track points towards a cluster, the cluster's energy is re- 
defined by subtracting the energy which is expected to be deposited in the calorimeter by the 
track. If the energy of the cluster is smaller than this expected energy, the cluster is not used. 
In this way, the energies of the clusters are primarily associated with neutral particles. 

Each accepted track and cluster is considered to be a particle. Tracks are assigned the pion 
mass. Clusters are assigned zero mass since they originate mostly from photons. 

To eliminate residual background and events in which a significant number of particles is 
lost near the beam direction, the number of accepted charged tracks in an event is required to 
be at least five and the thrust axis of the event, calculated using the particles, is required to 
satisfy | cos(# t hrast)| <0.90, where thrust is the angle between the thrust and beam axes. The 
number of events which passes these cuts is 2 407000. The residual background to this sample 
from all sources is estimated to be less than 1% [14] and is neglected. 



3 QCD models 

To establish the sensitivity of our analysis to processes with color reconnection, we generate 
events using Monte Carlo simulations of perturbative QCD and the hadronization process, both 
with and without the effects of reconnection. 

The models without reconnection in our study are the Jetset [16], Herwig [17, 18] and Ari- 
adne [19] Monte Carlo programs, versions 7.4, 6.2 and 4.11 respectively. Jetset and Herwig are 
based on parton showers with branchings described by Altarelli-Parisi splitting functions [20], 
followed by string hadronization [6] for Jetset and cluster hadronization [21] for Herwig. Ari- 
adne employs the dipole cascade model [22] to generate a parton shower, followed by string 
hadronization. The principal parameters of the models were tuned to yield an optimized de- 
scription of the global properties of hadronic Z° events and are documented in [23] for Jetset 
and in Tables 1 and 2 for Herwig and Ariadne. All three models provide a good description of 
the main features of e + e~ hadronic annihilation events, including the properties of identified 
gluon jets, see for example [11]. 

The models in our study which incorporate color reconnection are the model of Lonnblad [25] 
implemented in the Ariadne Monte Carlo 3 , the color reconnection model [18] in the Herwig 
Monte Carlo, and a model introduced by Rathsman [27]. We refer to these as the Ariadne-CR, 
Herwig-CR, and Rathsman-CR models, respectively. The Ariadne-CR model is an extension 
of the model of Gustafson and Hakkinen [28]. The Rathsman-CR model is implemented in 
the Pythia Monte Carlo [16], version 5.7. For e + e~ annihilations in the absence of initial-state 

3 There are three variants of the color reconnection model in Ariadne, corresponding to settings of the 
parameter MSTA(35) = 1, 2 or 3; for hard processes involving a single color singlet system, such as Z° decays, 
all three variants are identical; note that the parameter PARA(28) should be set to zero if the MSTA(35)=2 
option is used in Z° decays [26]. 
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Parameter 


Monte Carlo 
name 


Default 
value 


Optimized 
value 


Aqcd (GeV) 
Gluon mass (GeV/c 2 ) 
Maximum cluster mass parameter (GeV/c 2 ) 
Maximum cluster mass parameter 
Cluster spectrum parameter, udsc 

Cluster spectrum parameter, b 
Gaussian smearing parameter, udsc 
Decuplet baryon weight 


QCDLAM 
RMASS(13) 
CLMAX 
CLPOW 
PSPLT(l) 
PSPLT(2) 
CLSMR(l) 
DECWT 


0.18 

0.75 

3.35 

2.00 

1.00 

1.00 

0.0 

1.0 


0.18 ±0.01 
0.75 ±0.05 
3.35 ±0.05 
2.0 ±0.2 
1.00 ±0.05 

o ^ +0 - 07 

u -< 3 < 3 -0.03 

40+ - 20 
0.7 ±0.1 



Table 1: OPAL parameter set for Herwig, version 6.2. The method used to tune the parameters 
is presented in [23]. Parameters not listed were left at their default values. The uncertainties 
represent ±1 standard deviation limits obtained from the x 2 contours. The x 2 contours were 
defined by varying the parameters one at a time from their tuned values. 



Parameter 


Monte Carlo 
name 


Default 
value 


Optimized 
value 


Aqcd (GeV) 

PT,min. (GeV/c) 

b (GeV- 2 ) 
V(qq)/V(q) 
[V(us)/V(ud)]/[V(s)/V(d)] 

?W)/3?K) 


PARA(l) 
PARA (3) 
PARJ(42) 
PARJ(l) 
PARJ(3) 
PAR J (4) 


0.22 
0.60 
0.58 
0.10 
0.40 
0.05 


0.215 ±0.002 
0.70 ±0.05 
0.63 ±0.01 
0.130 ±0.003 
0.600 ±0.016 

U.U^U_ 003 


Extra Baryon suppression 
(MSTJ(12)=3) 


PARJ(19) 


1.00 


0.52 



Table 2: OPAL parameter set for Ariadne, version 4.11. The method used to tune the param- 
eters is presented in [23]. Parameters not listed were left at their default values. The extra 
baryon suppression factor PARJ(19), enabled by setting MSTJ(12)=3, was taken from [24]. 
The uncertainties have the same meaning as in Table 1. 
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photon radiation, Pythia is equivalent to Jetset. Thus, the Rathsman-CR model is effectively a 
version of Jetset which contains color reconnection. We note the Pythia Monte Carlo contains 
its own color reconnection model, based on the work of Khoze and Sjostrand [29]. We do not 
include this model in our study because it is not implemented for Z° decays. The Rathsman- 
CR model has been found to provide a good description of rapidity gap measurements in both 
electron-proton and proton-antiproton collisions [30]. 

The parameters we use for the Ariadne-CR model are the same as those given in Table 2 for 
Ariadne except for the parameter PARJ(42) which was adjusted from 0.63 to 0.55 GeV~ 2 so 
that the model describes the measured value of mean charged particle multiplicity in inclusive 
Z° decays, (n c h.), see Sect. 4. Analogously, the parameters of the Herwig-CR model are the same 
as those used for Herwig (see Table 1) except CLMAX was adjusted from 3.35 to 3.75 GeV/c 2 
and RMASS(13) from 0.75 to 0.793 GeV/c 2 to describe (n ch ). For our implementation of the 
Rathsman-CR model, we use the parameter set given for Jetset in [23]. 

Besides the Jetset parameters, the Rathsman-CR model employs a parameter, denoted Rq, 
which is an overall suppression factor for color reconnection. The value of Rq is not arbitrary 
but reflects the 1/Nq suppression of reconnected events compared to planar events mentioned 
in the Introduction. For this parameter, we use R = 0.1 as suggested in [27]. The analogous 
parameter in the Herwig-CR model, PRECO, is maintained at its default value of 1/Nq = 1/9. 
For the Ariadne-CR model, the corresponding parameter, PARA(26), stipulates the number of 
distinct dipole color states. We use the default value for this parameter, PARA(26) =9, which 
again corresponds to Nq = 3 and the l/N^ suppression of reconnected processes. 

Our implementations of the Ariadne-CR, Herwig-CR and Rathsman-CR models provide 
descriptions of the global features of e + e~ data which are essentially equivalent to those of the 
corresponding models without reconnection. This is discussed in Sect. 4 below. 

The Monte Carlo events are examined at two levels: the "detector level" and the "hadron 
level." The detector level includes initial-state photon radiation, simulation of the OPAL 
detector [31], and the analysis procedures described in Sect. 2. The hadron level does not 
include these effects and utilizes all charged and neutral particles with lifetimes greater than 
3 x 10~ 10 s, which are treated as stable. Samples of 6 million Ariadne, Herwig and Jetset events, 
and 3 million Ariadne-CR, Herwig-CR and Rathsman-CR events, were processed through the 
detector simulation and used as the detector level samples in this study. The hadron level 
samples are based on 10 million Monte Carlo events for each model. 



4 Model predictions for inclusive Z u decays 

The Ariadne-CR, Rathsman-CR and Herwig-CR models yield descriptions of standard measures 
of properties in inclusive Z° data which are essentially equivalent to those provided by Ariadne, 
Jetset, and Herwig, respectively, as stated above. Thus, color reconnection as implemented in 
these models has only a small effect on the global features of inclusive e + e~ events. To illustrate 
these points, we measured the following distributions using the inclusive Z° sample discussed 
in Sect. 2: 
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1. Sphericity, S [32]; 



2. Aplanarity, A [32]; 

3. the negative logarithm of the jet resolution scale for which an event changes from being 
classified as a three-jet event to a four-jet event, using the Durham jet finder [33], — ln(y 34 ); 

4. charged particle rapidity with respect to the thrust axis, yr- 

Note there are correlations between these variables and between different bins of some of the 
distributions. Using a sample of Jetset events at the hadron level, the correlation coefficient 
between S and A was found to be 0.66, between S and — ln(y 34 ) —0.61, and between A and 
— ln(y 34 ) —0.66. Similar results were found using the other models. The yr distribution contains 
one entry per particle, in contrast to the other distributions which contain one entry per event. 
Therefore, the yx distribution was not included in this correlation study. Taken together, the 
four distributions are sensitive to the momentum structure of an event both in and out of the 
three-jet event plane, to four-jet event structure, and to particle multiplicity. They therefore 
provide a relatively complete and relatively uncorrelated set of distributions with which to 
assess the global features of e + e~ events. 

The distributions were corrected to the hadron level using bin-by-bin factors. The method 
of bin-by-bin corrections is described in [32]. Ariadne was used to determine the correction 
factors. Ariadne was chosen because it was found to provide a better description of the data 
at the detector level than Jetset or Herwig. The typical size of the corrections is 10%. As 
systematic uncertainties, we considered the following. 

• The other models - Jetset, Herwig, Ariadne-CR, Herwig-CR and Rathsman-CR - were 
used to determine the correction factors, rather than Ariadne. 

• Charged tracks alone were used for the data and Monte Carlo samples with detector 
simulation, rather than charged tracks plus electromagnetic clusters. 

• The particle selection was further varied, first by restricting charged tracks and elec- 
tromagnetic clusters to the central region of the detector, |cos#| < 0.70, rather than 
| cos^l < 0.96 for the charged tracks and | cos#| < 0.98 for the clusters, and second by 
increasing the minimum transverse momentum of charged tracks with respect to the beam 
axis from 0.05 GeV/c to 0.15 GeV/c. 

The differences between the standard results and those found using each of these conditions were 
used to define symmetric systematic uncertainties. For the first item, the largest of the described 
differences with respect to the standard result was assigned as the systematic uncertainty, and 
similarly for the third item. The systematic uncertainties were added in quadrature to define 
the total systematic uncertainties. The systematic uncertainty evaluated for each bin was 
averaged with the results from its two neighbors to reduce the effect of bin-to-bin fluctuations. 
The single neighbor was used for bins at the ends of the distributions. 

The largest contribution to the systematic uncertainty of the S, — ln(y 34 ) and yx distribu- 
tions arose from using the Herwig-CR model to correct the data. For A, the largest systematic 
effect was from using Jetset to correct the data. 
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Model 


S 


A 


- In (1/34) 


Vt 


Total 


(rich.) 


(Number of bins) 


(19) 


(15) 


(26) 


(21) 


(81) 




AriflHnp 

ill J.Q I. i.±±v3 


1 7 


6 8 


10 7 


17 7 


36 9 


21 06 


Ariadne- CR 


6.2 


5.9 


4.3 


16.0 


32.4 


21.09 


Jet set 


18.4 


91.4 


64.1 


26.8 


200.7 


21.09 


Rathsman-CR 


18.5 


103.6 


74.7 


46.7 


243.5 


20.80 


Herwig 


19.3 


27.0 


42.5 


39.1 


127.9 


21.14 


Herwig-CR 


in c 

lu.o 




on p. 

0U.0 


y4.o 


ioi.o 


01 nfi 
zl.Uo 


Re-tuned Ariadne-CR 


498.1 


548.7 


1001.3 


971.2 


3019.3 


21.12 


(pr,min. = 4.7 GeV/c, 6 = 0.17 GeV" 2 ) 














Re-tuned Rathsman-CR 


106.8 


294.2 


429.7 


287.0 


1117.7 


21.16 


(Qo = 5.5 GeV/c 2 , 6 = 0.27 GeV" 2 ) 















Table 3: x 2 values between the data and models for the distributions shown in Figs. 2-5, 
calculated using the full experimental uncertainties including systematic terms. The number of 
bins in each distribution is given in parentheses in the second row. The last two rows give the 
results for re-tuned versions of the Ariadne-CR and Rathsman-CR models, see Sect. 7. 4. The 
models' predictions for the mean charged particle multiplicity in inclusive Z° decays, (n c h.), are 
listed in the last column. 



The corrected measurements of S, A, — ln(y 34 ) and yx are presented in Figs. 2-5. These data 
are consistent with our previously published results [32]. The data are shown in comparison to 
the predictions of the models at the hadron level. The model predictions are generally seen to be 
similar to each other and in agreement with the experiment. Parts (b) and (c) of Figs. 2-5 show 
the deviations of the Monte Carlo predictions from the data in units of the total experimental 
uncertainties "odata," with statistical and systematic terms added in quadrature. The statistical 
uncertainties are negligible compared to the systematic uncertainties. The curves labelled 
"Re-tuned Rathsman-CR" and "Re-tuned Ariadne-CR" in parts (b) and (c) are discussed in 
Sect. 7.4. 

We calculated the \ 2 values between the hadron level predictions of the models and the 
corrected data. The x 2 values were determined using the total experimental uncertainties, with 
no accounting for correlations between the different bins or distributions. The x 2 results are 
listed in Table 3. These x 2 values are intended to be used only as a relative measure of the 
description of the data by the models. Since the uncertainties are dominated by systematics and 
correlations are not considered, these x 2 values cannot be used to determine confidence levels 
assuming the uncertainties are distributed according to a normal distribution. In particular, a 
good description does not imply that a model's x 2 should approximately equal the number of 
data bins. 

From Table 3, it is seen that the x 2 results for Ariadne are much smaller than for Jetset or 
Herwig. The reason for this is partly that the detector level distributions are better described 
by Ariadne, as stated above, and partly that Ariadne is used to determine the correction factors. 
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It is unavoidable that the correction procedure introduces a bias towards the model used to 
perform the corrections, as discussed for example in [32]. These biases - although small - can 
have a significant effect on the x 2 values because of the small experimental uncertainties. For 
this reason, it is only meaningful to compare our x 2 results within the context of a specific 
parton shower and hadronization scheme, e.g. Ariadne with Ariadne-CR but not Ariadne with 
Jet set. 

The total x 2 f°r the Ariadne-CR model is seen to be about the same as for Ariadne (in 
fact it is a little smaller). For the Herwig-CR model, the total x 2 is about 20% larger than for 
Herwig. From Table 3, it is seen that this difference arises entirely from a single distribution, y T) 
however. Similarly, the total x 2 f° r the Rathsman-CR model is about 20% larger than for Jetset, 
with the largest contribution to the increase from y T . 

The last column in Table 3 lists the predictions of the models for the mean value of charged 
particle multiplicity, (n c h.). The statistical uncertainties are negligible. The results of all 
models agree with the LEP-averaged result for Z° decays, (n ch ) = 21.15 ± 0.29 [34], to within 
the uncertainties. 

Thus the three models with color reconnection yield overall descriptions of the global proper- 
ties of hadronic Z° events which are essentially equivalent to those of the corresponding models 
without reconnection. For the Ariadne-CR and Rathsman-CR models, this agrees with the 
observations in [25] and [27], respectively. 



5 Gluon jet selection 

To define jets, we use the Durham jet finder [33]. The resolution scale, y cnt , is adjusted sepa- 
rately for each event so that exactly three jets are reconstructed. The jets are assigned energies 
using the technique of calculated energies with massive kinematics, see for example [11]. This 
method relies primarily on the angles between jets and the assumption of energy-momentum 
conservation. Jet energies determined in this manner are more accurate than visible jet ener- 
gies, with the latter defined by a sum over the reconstructed energies of the particles assigned 
to the jet. 

To identify which of the three jets is the gluon jet, we reconstruct displaced secondary 
vertices in the quark (q or q) jets and thereby anti-tag the gluon jet. Displaced secondary 
vertices are associated with heavy quark decay, especially that of the b quark. At LEP, b quarks 
are produced almost exclusively at the electroweak vertex 4 : thus a jet containing a b hadron 
is almost always a quark jet. To reconstruct secondary vertices in jets, we employ the method 
described in [36]. Briefly, charged tracks are selected for the secondary vertex reconstruction 
procedure if they are assigned to the jet by the jet finder, have coordinate information from 
at least one of the two silicon detector layers, a momentum larger than 0.5 GeV/c, and a 
distance of closest approach to the primary event vertex [36] less than 0.3 cm. In addition, the 
uncertainty on the distance of closest approach must be less than 0.1 cm. A secondary vertex 
is fitted using the so-called "tear down" method [36] and is required to contain at least three 

4 About 22% of hadronic Z° events contain a bb quark pair from the electroweak decay of the Z°, compared 
to only about 0.3% [35] with a bb pair from gluon splitting. 
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Figure 2: (a) The Sphericity distribution for inclusive Z° events, in comparison to the predic- 
tions of models with and without color reconnection (CR). The data have been corrected for 
initial-state photon radiation and detector response. The statistical uncertainties are too small 
to be visible. The vertical lines attached to the data points (barely visible) show the total 
uncertainties, with statistical and systematic terms added in quadrature, (b) and (c) show 
the deviations of the Monte Carlo predictions from the data in units of the total experimental 
uncertainties, Odata- 
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Figure 3: (a) The Aplanarity distribution for inclusive Z° events, in comparison to the pre- 
dictions of models with and without color reconnection (CR). The data have been corrected 
for initial-state photon radiation and detector response. The uncertainties - both statistical 
and total - are too small to be visible, (b) and (c) show the deviations of the Monte Carlo 
predictions from the data in units of the total experimental uncertainties, Odata- 



13 



0.4 



^ 0.3 



o.i - 



13 







10 



a 



-10 
10 



~i — i — i — i — i — i — i — i — i — i — r 

$ OPAL 

Jetset7.4 

Rathsman-CR 

Ariadne 4.11 

— — Ariadne- CR 

Herwig 6.2 

Herwig-CR 



~i — i — i — i — i — i — i — i — r 



l I ■ ■ J JW>i 



o 




(a) 



10 



12 



Jetset 7.4 
Rathsman-CR 



(b) 




Re-tuned Rathsman-CR 
Herwig 6.2 



J I L 



J I L 



J I L 



J I L 



J L 



J I L 







8 



10 



12 



S3 



5 - 



I 



-5 - 



-10 



i — i — r 



i — i — r 

Ariadne 4.11 
Ariadne-CR 



i — i — r 



7^ 



i — i — r 



i — i — r 



i — i — r 



(c) 



• Re-tuned Ariadne-CR 
- Herwig-CR 



_l I I La I I l_ 



i i i 







8 



10 



12 



Figure 4: (a) The — ln(2/ 34 ) distribution for inclusive Z° events, in comparison to the predictions 
of models with and without color reconnection (CR). The data have been corrected for initial- 
state photon radiation and detector response. The uncertainties - both statistical and total - 
are too small to be visible, (b) and (c) show the deviations of the Monte Carlo predictions from 
the data in units of the total experimental uncertainties, o" data . 
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Figure 5: (a) The yT distribution for inclusive Z° events, in comparison to the predictions of 
models with and without color reconnection (CR). The data have been corrected for initial- 
state photon radiation and detector response. The uncertainties - both statistical and total - 
are too small to be visible, (b) and (c) show the deviations of the Monte Carlo predictions from 
the data in units of the total experimental uncertainties, Odata- 
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such tracks. For jets with such a secondary vertex, the signed decay length, L, is calculated 
with respect to the primary vertex, along with its error, a^. The sign of L is determined by 
summing the 3-momenta of the tracks fitted to the secondary vertex; L > if the secondary 
vertex is displaced from the primary vertex in the same hemisphere as this momentum sum, 
and L < otherwise. To be identified as a quark jet, a jet is required to have a successfully 
reconstructed secondary vertex with L/a L > 2.0 if it is the highest energy jet or L/a L > 5.0 if 
it is one of the two lower energy jets. We require the highest energy jet and exactly one of the 
two lower energy jets to be identified as quark jets. The other lower energy jet is tagged as a 
gluon jet. 

For each tagged gluon jet, we determine the scale, /tj et , given by 



where Ej et is the energy of the jet, with # min . the smaller of the angles between the gluon jet 
and the other two jets. Note that due to QCD coherence, the properties of a gluon jet in 
e + e~ annihilations depend on a transverse momentum-like quantity such as n- iCt and not the jet 
energy, see for example [37]. Kj et as defined in eq. (1) was shown to be an appropriate scale for 
gluon jets in [38]. 

The Kjet distribution of the tagged gluon jets is shown in Fig. 6. The data are presented in 
comparison to the predictions of the detector level QCD models introduced in Sect. 3. All the 
simulations are seen to provide a good description of the measured n- )Ct spectrum. 

To select hard, acollinear gluon jets, we require /Cj et > 7 GeV. Further, we require the energy 
of the gluon jets to be less than 35 GeV because the simulations predict the gluon jet purity 
(see below) drops sharply for higher energies. The jets are required to contain at least two 
particles. With these cuts, the number of selected gluon jets is 12 611. The energy of the jets 
varies from about 10 GeV up to the cutoff of 35 GeV, with an average and RMS of 21.7 GeV 
and 6.6 GeV, respectively. 

To evaluate the purity of the gluon jets, we use Monte Carlo samples at the detector level. 
We determine the directions of the primary quark and antiquark from the Z° decay after the 
parton shower has terminated. The detector level jet closest to the direction of an evolved 
primary quark or antiquark is considered to be a quark jet. The distinct jet closest to the 
evolved primary quark or antiquark not associated with this first jet is considered to be the 
other quark jet. The remaining jet is the gluon jet. The estimated gluon jet purity found using 
Jetset is approximately constant at 98% for jet energies from 10 to 25 GeV, then decreases to 
78% at 35 GeV. The overall purity is (94.6 ± 0.1 (stat.))%. Similar results are obtained using 
all other models except for Ariadne-CR. 5 Note the overall purity of the gluon jets decreases to 
86% after the requirement of a rapidity gap is imposed, see Sect. 6. The reason the purity is 
lower if a rapidity gap is required is because gluon jets have a larger mean multiplicity than 
quark jets [39], making it less likely a gap will occur in gluon jets compared to quark jets as 
the result of a fluctuation. By requiring the presence of a rapidity gap, the relative proportion 
of quark jets is therefore enhanced. 

5 For the Ariadne-CR model, the estimated purity is smaller, about 72%. Since this model does not describe 
our gluon jet measurements well (see Sect. 7), it is not clear if this estimate is reliable, however. Note the 
estimates of gluon jet purity are presented for informational purposes only. 




(1) 
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Figure 6: Distribution of the n- iCt scale of tagged gluon jets, see eq. (1). The distribution 
includes the effects of initial-state photon radiation and detector acceptance and resolution. The 
uncertainties are statistical only. The results are shown in comparison to the predictions of QCD 
Monte Carlo programs which include detector simulation and the same analysis procedures as 
are applied to the data. To define hard, acollinear gluon jets, the region /tj Ct > 7 GeV, to the 
right of the vertical dashed line, is selected. The hatched area shows the quark jet background 
evaluated using Jetset. 

6 Rapidity gap analysis 



To identify gluon jets with a rapidity gap, we examine the charged and neutral particles assigned 
to the selected gluon jets by the jet finder. The rapidities of the particles are determined with 
respect to the jet axis. The particles in the jet are ordered by their rapidity values. 

Models with color reconnection are expected to yield more events with a large rapidity gap 
than models without reconnection, as discussed in the Introduction. A large rapidity gap can 
correspond to a large value for the smallest particle rapidity in a jet, y m i n , or else to a large 
value for the maximum difference between the rapidities of adjacent rapidity-ordered particles, 
Ay max . These two types of rapidity gap conditions are illustrated schematically in Fig. 7. Note 
that the Durham jet finder occasionally assigns particles to a jet even if the angle between the 
particle and jet axis is greater than 90°. This explains the negative rapidity values illustrated 
for some particles in Fig. 7b. 

The measured distribution of 2/ min is presented in Fig. 8a. The data are shown in comparison 
to the predictions of the models at the detector level. To emphasize the difference between 
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(a) 1 dn 



(b) 




max 



Figure 7: Schematic illustration of the distribution of particle rapidities for gluon jets with a 
rapidity gap as defined in this study: (a) for the y m i n sample (see text), and (b) for the Ay max 
sample. The leading parts of the gluon jets are defined by charged and neutral particles with 
rapidities beyond the gap, as indicated in the figure. 
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models with and without color reconnection, we form the following ratio: 

x /(2/min)cR - fiVmrn) noCR / c\\ 

"j/min — 77~ \ \ Z ) 

/ U/minJno CR 

with /(y m in)cR the prediction of a model with color reconnection for a bin of the y min distribution 
in Fig. 8a and /(y m in)noCR the prediction of the corresponding model without reconnection. 
The results for S yinin are shown in Fig. 8b. For the Rathsman-CR model, a significant excess of 
events is observed relative to Jetset for y m i n values larger than about 1.4, and similarly for the 
Ariadne-CR model relative to Ariadne. The Herwig-CR model exhibits a similar excess with 
respect to Herwig, although with less significance. Based on these results, we choose y min > 1.4 
to select a sample of gluon jets with a rapidity gap, see the dashed vertical line in Fig. 8b. In 
the following, we refer to this as the "ymin" sample. 

For gluon jets with y mm <l-4, we measure Ay max . The resulting distribution is shown in 
Fig. 9a. In analogy to eq. (2), we form the fractional difference S^y^- The distribution of 
$A?/ ma x is shown in Fig. 9b. A significant excess of events is observed for the Ariadne-CR and 
Rathsman-CR models, relative to Ariadne and Jetset, for Ay max larger than about 1.3. We 
therefore choose Ay max > 1.3 to select an additional sample of gluon jets with a rapidity gap, 
see the dashed vertical line in Fig. 9b. In the following, we refer to this as the "Ay max " sample. 
For the Herwig-CR model, there is not a clear excess of events relative to Herwig for any 
Ay max value, suggesting this distribution is not sensitive to color reconnection as implemented 
in Herwig. In the following, we therefore test the Herwig-CR model using the y m i n sample only, 
not the standard data set defined by the y min and Ay max samples taken together. 

In total, 655 gluon jets with a rapidity gap are selected, 496 in the y m m sample and 159 
in the Ay max sample. The purity of the gluon jets, evaluated using the method described in 
Sect. 5, is approximately 94% for gluon jet energies between 10 and 25 GeV, then drops to 
about 50% at 35 GeV. The overall purity is (85.7 ± 1.0 (stat.))%. Our subsequent study is 
based on the leading part of these jets, defined by charged and neutral particles with y > y min 
for events 6 in the y m i n sample and by particles with rapidities beyond the gap A|/ max for events 
in the Ay max sample, see Fig. 7. 



7 Color reconnection study 

To remain as sensitive as possible to color reconnection, we first compare the Monte Carlo 
distributions to the data at the detector level. Following this, we correct the measurements for 
the effects of initial-state radiation, detector acceptance and resolution, and gluon jet impurity, 
and compare the predictions of the models to the data at the hadron level. The hadron level 
study allows us to more readily assess the effect of adjusting Monte Carlo parameters, see 
Sect. 7.4. 

The distributions presented in this section are normalized to the total number of selected 
gluon jets discussed in Sect. 5, i.e. to the number of gluon jets before the rapidity gap require- 
ment. The reason for this is to remain sensitive to the rate at which gluon jets with a rapidity 
gap occur, e.g. to the production rate of events like Fig. lb. 

6 For events in this class, this is therefore the entire jet, see Fig. 7a. 
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Figure 8: (a) Distribution of y min in the tagged gluon jets. The distribution includes the effects 
of initial-state photon radiation and detector acceptance and resolution. The uncertainties 
are statistical only. The results are shown in comparison to the predictions of QCD Monte 
Carlo programs which include detector simulation and the same analysis procedures as are 
applied to the data. The hatched area shows the quark jet background evaluated using Jetset. 
(b) Fractional difference between the results of a Monte Carlo program with color reconnection 
and the corresponding model without reconnection. To define gluon jets with a rapidity gap, 
the region 7/ min > 1.4, to the right of the vertical dashed line, is selected. 
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Figure 9: (a) Distribution of Ay max for tagged gluon jets with y min <1.4. The distribution 
includes the effects of initial-state photon radiation and detector acceptance and resolution. The 
uncertainties are statistical only. The results are shown in comparison to the predictions of QCD 
Monte Carlo programs which include detector simulation and the same analysis procedures as 
are applied to the data. The hatched area shows the quark jet background evaluated using 
Jetset. (b) Fractional difference between the results of a Monte Carlo program with color 
reconnection and the corresponding model without reconnection. To define gluon jets with a 
rapidity gap, the region Ay max > 1.3, to the right of the vertical dashed line, is selected. 
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7.1 Detector level distributions 



The charged particle multiplicity distribution of the leading part of the gluon jets, fading) 
is shown in Fig. 10a. The results are shown in comparison to the predictions of the Jet- 
set and Rathsman-CR models. Fig. 10b shows the same data compared to Ariadne and 
Ariadne-CR. The most striking feature of these results is the large excess of entries predicted by 
the Ariadne-CR and Rathsman-CR models at nfg ading = 2 and 4 compared to the corresponding 
models without color reconnection. Using Monte Carlo information, we verified these excesses 
are a consequence of events like Fig. lb, present in the CR models but not in the models with- 
out CR. The isolated, electrically neutral gluonic system in the leading part of the gluon jets in 
these events decays into an even number of charged particles, yielding the spikes at nf^ ading = 2 
and 4. 

The data are generally well described by Jetset (Fig. 10a), except for the bins with nfg ading = 
1, 2 and 4 where the data exceed the predictions by more than one standard deviation of the 
statistical uncertainties. The description by Ariadne (Fig. 10b) is considerably worse in that 
the data lie well above the Ariadne results for most of the range between n^ ding = 2 and 6. 
Nonetheless, Jetset and Ariadne provide a much better overall description of the data than 
the corresponding models with reconnection. In particular, there is not a significant "spiking 
effect" in the data at even values of multiplicity as predicted by these two CR models. We 
conclude that color reconnection as implemented by the Rathsman-CR and Ariadne-CR models 
is strongly disfavored, at least using their standard parameters given in Sect. 3. 

The ^i^ing distribution obtained using the y min selection (see Sect. 6) is presented in Fig. 11. 
The data are shown in comparison to the corresponding results of the Herwig and Herwig-CR 
models. We use the y min selection to test Herwig-CR, and not the standard selection defined by 
the combined y min and Ay max samples, because the latter is not sensitive to differences between 
the Herwig and Herwig-CR models as discussed in Sect. 6. For purposes of comparison, the 
prediction of Herwig using the standard selection is shown in Fig. 10a, however. 

From Fig. 11, the Herwig-CR model is seen to predict a systematic excess of entries relative 
to the corresponding model without CR for multiplicities between about 2 and 5. The overall 
description of the nf^ adi distribution by the Herwig-CR model is nonetheless reasonable, at 
least in comparison to the predictions of Jetset and Ariadne in Fig. 10. The best overall 
description of the n^ ding distribution is provided by Herwig. 

We next sum the charges of the particles in the leading part of the gluon jets to find the total 
leading electric charge, Qieading- This type of distribution was suggested in [9]. The distribution 
of Qieading is shown in Fig. 12a. The Rathsman-CR and Ariadne-CR models are seen to predict 
a large excess of events with Qieading = compared to the data or models without reconnection, 
due to the presence of electrically neutral isolated gluonic systems at large rapidities as discussed 
above. The Jetset and Ariadne predictions for the rate of gluon jets with Qieading — are about 
20% too low. For purposes of comparison, the prediction of Herwig is shown in Fig. 12a. Herwig 
is seen to describe the data well. 

Our measurement of the rate of gluon jets with (pleading = therefore lies between the predic- 
tions of the Jetset and Rathsman-CR models, and similarly between the predictions of Ariadne 
and Ariadne-CR. In this respect, the data appear to be consistent with the presence of a fi- 
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Figure 10: Distribution of n£r ading in the leading part of gluon jets, based on our standard 
selection. "N" represents the total number of selected gluon jets and "N gap " the number of 
gluon jets with a rapidity gap. The distribution includes the effects of initial-state photon 
radiation and detector acceptance and resolution. The uncertainties are statistical only. The 
results are shown in comparison to the predictions of QCD Monte Carlo programs which include 
detector simulation and the same analysis procedures as are applied to the data: (a) the Jetset, 
Rathsman-CR and Herwig models, and (b) the Ariadne and Ariadne-CR models. The hatched 
area shows the quark jet background evaluated using Jetset in part (a) and Ariadne in part (b). 
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Figure 11: Distribution of n^ ading in the leading part of gluon jets, based on the y min selection 
(see Sect. 6). "N" represents the total number of selected gluon jets and "N gap " the number 
of gluon jets with a rapidity gap. The distributions include the effects of initial-state photon 
radiation and detector acceptance and resolution. The results are shown in comparison to the 
predictions of the Herwig and Herwig-CR models. The uncertainties are statistical only. The 
hatched area shows the quark jet background evaluated using Herwig. 



nite amount of color reconnection, at least as predicted by these two CR models, albeit at a 
significantly smaller level than predicted by the default CR settings of the models. The most 
unambiguous signal for color reconnection in our study is the spiking effect at even values 
of Jading seen m Fig- 10, however. The data do not provide clear evidence for these spikes. 
Furthermore, the Herwig model without CR describes the (pleading distribution well, as seen 
from Fig. 12a. Therefore, the discrepancies of Jetset and Ariadne with the data in Fig. 12a 
do not provide unambiguous evidence for reconnection effects, but instead are consistent with 
other inadequacies in the simulations, not related to CR. The same statement holds for the 
discrepancies of Jetset and Ariadne with the data in Fig. 10. 

The distribution of Qi ea ding obtained using the y min selection is presented in Fig. 12b. The 
data are shown in comparison to the corresponding results from Herwig and Herwig-CR. Herwig 
describes the data well, similar to Fig. 12a. The predictions of the Herwig-CR model are seen 
to lie somewhat above the data, especially for fading = and 1. 

As a systematic check, we repeated the analysis presented above using different choices for 
the scale of gluon jets, /Cj et , see eq. (1). Specifically we examined the results for 4 < Kj et < 7 GeV 
and Kj e t < 4 GeV, rather than /tj Ct > 7 GeV as in our standard analysis. Note the definition of 
a gluon jet becomes ambiguous for small Kj Ct values. We find that the spikes at even values of 
n kadi predicted by the Rathsman-CR and Ariadne-CR models become much less prominent 
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Figure 12: Distribution of Qieadmg in the leading part of gluon jets in comparison to the pre- 
dictions of QCD Model Carlo programs: (a) using the standard selection; (b) using the y m i n 
selection. "N" represents the total number of selected gluon jets and "N gap " the number of 
gluon jets with a rapidity gap. For both parts (a) and (b), the distributions include the effects 
of initial-state photon radiation and detector acceptance and resolution. The uncertainties are 
statistical only. The hatched area shows the quark jet background evaluated using Herwig. 
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for the smaller «j et scales, especially the spike at nf eading = 4, i.e. the selections with softer 
or more collinear gluon jets are less sensitive to color reconnection. This justifies the choice 
Kjet > 7 GeV for our standard analysis. To the extent that a CR signal is still visible using the 
smaller /tj Ct ranges, we find that the values of y m i n and A?/ max above which the predictions of the 
CR models exhibit deviations from the non-CR models are similar to those shown in Figs. 8b 
and 9b. 

As an additional check, we repeated the analysis described in Sects. 5 and 6 except using 
energy ordering to identify gluon jets rather than secondary vertex reconstruction. In the 
energy ordering method, the jet with the smallest calculated energy in three-jet qqg events is 
assumed to be the gluon jet. The purity of gluon jets identified using this technique is much 
lower than found using secondary vertices, especially for the high energy jets most sensitive to 
color reconnection. The gluon jet purity found using energy ordering is 64%, compared to 95% 
for our standard analysis. To increase the purity, we therefore required £j et < 15 GeV, rather 
than £j et < 35 GeV as in the standard analysis. This method yields about 94 000 tagged gluon 
jets. The mean gluon jet energy is 12.9 GeV and the estimated purity 81%. After imposing the 
rapidity gap requirements of Sect. 6, we obtain 6604 gluon jets with an estimated purity of 56%. 
The results we obtain from this check are consistent with our observations presented above. 
In particular, the results for the Qi ea ding distribution are qualitatively similar to those shown 
in Fig. 12. We note, however, that the spike at n^ ading = 4 predicted by the Rathsman-CR 
and Ariadne-CR models in Fig. 10 is not visible in the corresponding Monte Carlo predictions 
based on energy ordering, because of their softer energy scales (i.e. this is similar to the check 
employing smaller /tj Ct values, mentioned above). Therefore the selection using energy ordering 
is not as sensitive to color reconnection as our standard selection. 

The results of Figs. 10-12 demonstrate the sensitivity of our study to processes with color 
reconnection. We discuss the effect of parameter variation on the predictions of the Rathsman- 
CR and Ariadne-CR models in Sect. 7.4. 



7.2 Correction procedure 

As the next step in our study, we correct the data in Figs. 10-12 to the hadron level. The 
correction procedure employs an unfolding matrix. The matrix is constructed using detector 
level Monte Carlo events. The events are subjected to the detector level requirements of Sects. 2 
and 5. In addition, the events are required to exhibit a rapidity gap, defined by the conditions 
of Sect. 6, at both the detector and hadron levels. The matrices relate the values of Tiding 
and Qicading at the detector level to the corresponding values before the same event is processed 
by the detector simulation. Thus the matrices correct the data to the hadron level with the 
exception that initial-state radiation and the experimental event acceptance are included. In 
a second step, the data are corrected for event acceptance, initial-state radiation, and gluon 
jet impurity using bin-by-bin factors. The matrices and bin-by-bin factors are determined 
using Herwig because the data in Figs. 10-12 are best described by that model. Statistical 
uncertainties are evaluated for the corrected data using propagation of errors, including the 
statistical uncertainties of the correction factors. 

Because of finite acceptance, especially for soft particles, significantly more events satisfy the 
rapidity gap requirements at the detector level than at the hadron level. As a consequence, the 
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overall corrections are fairly large, of the order of 40%. To verify the reliability of the correction 
procedure, we therefore performed the following test. We treated our sample of Jetset events 
at the detector level as "data," using the Herwig derived corrections to unfold them. The 
corrected Jetset distributions were found to agree with the corresponding Jetset distributions 
generated at the hadron level to within the statistical uncertainties. This demonstrates that 
our correction procedure does not introduce a significant bias. 

To evaluate systematic uncertainties for the corrected data, we repeated the analysis using 
the three systematic variations listed in Sect. 4, with one exception: to determine the systematic 
uncertainty related to the model dependence of the correction factors, we repeated the analysis 
using the Jetset, Ariadne and Herwig-CR models only. We did not include the Rathsman-CR 
or Ariadne-CR model because of their poor description of the data, see Figs. 10 and 12a. In 
addition, we made the following change to the standard analysis to assess the effect of altering 
the criteria used to identify gluon jets. 

• To identify the lower energy quark jets, we required the decay length to satisfy L/a L > 3 
rather than L/ctl > 5; this resulted in 1002 tagged gluon jets which satisfied the rapidity 
gap requirements, with an estimated purity of 76%. 

The systematic uncertainties were treated as described in Sect. 4, i.e. the full differences of the 
results of the systematic checks with respect to the standard analysis defined the systematic 
uncertainty for each term, and the individual terms were added in quadrature to define the 
total systematic uncertainties. 

The largest contributions to the total systematic uncertainties arose from using Ariadne to 
determine the correction factors, followed by the requirement L/a L > 3 to identify the lower 
energy quark jets. 



7.3 Hadron level distributions 

The corrected distributions of ^ ading and (pleading are presented in Fig. 13. These results are 
based on our standard selection, i.e. the y m i n and Ay max samples added together. The data 
are shown in comparison to the hadron level predictions of the Jetset, Rathsman-CR, Ariadne 
and Ariadne-CR models. For purposes of comparison, the predictions of Herwig are shown 
as well. The qualitative features of the predictions are seen to be similar to those of the 
corresponding detector level distributions in Figs. 10 and 12a. In particular, the Ariadne-CR 
and Rathsman-CR models exhibit a large excess of entries at n^ ading = 2 and 4 in Fig. 13a, 
corresponding to the <3ieading = bin in Fig. 13b, analogous to the results of Sect. 7.1. From 
Fig. 13a it is also seen that the Rathsman-CR model predicts a spike at nf^ ading = 6. This latter 
feature was not apparent in the detector level distribution of Fig. 10a because of finite detector 
resolution. 

The corresponding results based on the y min sample are presented in Fig. 14. The data are 
shown in comparison to the predictions of the Herwig and Herwig-CR models. From Fig. 14a 
it is seen that the Herwig-CR model predicts a significant excess of events relative to Herwig 
for n^ adi =2, 4 and 6, analogous to the results of the Rathsman-CR and Ariadne-CR models 
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in Fig. 13a. This suggests that the production of events like Fig. lb is a general feature of color 
reconnection. The spike in the prediction of the Herwig-CR model at n^ ading =4 in Fig. 14a 
probably explains the general excess of the Herwig-CR results above Herwig for multiplicities 
between nf^ ading =3 and 5 in Fig. 11. 

7.4 Effect of parameter variation on the model predictions 

We next study the effect of parameter variation on the predictions of the Rathsman-CR and 
Ariadne-CR models, to determine if they can be tuned to describe the gluon jet data of Fig. 13 
without adversely affecting their descriptions of the inclusive Z° decay measurements presented 
in Sect. 4. 

To begin, we define AQ\^^g ta to be the difference between the Monte Carlo prediction and 
experimental result for the Q leading — bin in Fig. 13b. We then vary the principal parameters 
of the models one at a time, with the other parameters at their standard values, to see if it 
is possible to reduce AQSh^* 11 to zero or near-zero, i.e. to obtain agreement of the model's 
prediction with this measurement. We note that if AQ^^[^ ta is near-zero, the predictions of 
the model for the fading distribution in Fig. 13a will also be in general agreement with the 
data since the events which yield the excess of entries in Fig. 13b are the same as those which 
yield the excess in Fig. 13a. 

The Rathsman-CR model: For the Rathsman-CR model, the following parameters were 
varied: 

• Aqcd, the QCD scale parameter, given by PARJ(81); 

• Qo, the minimum mass value to which partons evolve, given by PARJ(82); 

• a and b, which control the longitudinal momentum spectrum of hadrons relative to the 
string direction in the Lund model of hadronization, given by PARJ(41) and PARJ(42); 

• cr q , which controls the transverse momentum spectrum of hadrons, given by PARJ(21). 

The PARJ references are the names of the parameters in the Pythia Monte Carlo. These five 
parameters are the most important ones controlling the multiplicity and momentum distribu- 
tions of hadrons in the model. Note we do not include the color reconnection suppression factor 
i? mentioned in Sect. 3 in the above list. It is a trivial result that the Rathsman-CR model 
will describe the data as well as Jetset for R — > since the two models are identical in this 
limit. Varying R to reproduce the experimental result and corresponding uncertainty for the 
Qicading = bin in Fig. 13b yields R = 0.0085 ± 0.0075 (stat.) ± 0.0087 (syst.), consistent with 
R = 0. 

The results for AQ^-data 

are shown in Fig. 15. The standard values of the parameters 
are indicated by solid dots. The uncertainties attributed to the parameter values in [23], 
beyond which the description of inclusive Z° measurements is significantly degraded if the 
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Figure 13: Distributions of (a) n* ading and (b) Qieading in the leading part of gluon jets, based 
on our standard selection. "N" represents the total number of selected gluon jets and "N gap " 
the number of gluon jets with a rapidity gap. The data have been corrected for initial-state 
photon radiation, gluon jet impurity, and detector response. The horizontal bars indicate the 
statistical uncertainties. The vertical lines show the total uncertainties, with statistical and 
systematic terms added in quadrature. The results are shown in comparison to the predictions 
of QCD Monte Carlo models. 
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Figure 14: Distributions of (a) n^ ading and (b) Qicading in the leading part of gluon jets, based 
on the 2/ m i n selection. "N" represents the total number of selected gluon jets and "N gap " the 
number of gluon jets with a rapidity gap. The data have been corrected for initial-state photon 
radiation, gluon jet impurity, and detector response. The horizontal bars indicate the statistical 
uncertainties. The vertical lines show the total uncertainties, with statistical and systematic 
terms added in quadrature. The results are shown in comparison to the predictions of the 
Herwig and Herwig-CR Monte Carlo models. 
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other parameters remain at their standard values, are indicated by the horizontal error ranges. 
Note that an uncertainty is not evaluated for the a parameter in [23] and that the uncertainties 
attributed to Aqcd and b are too small to be visible. The width of the shaded bands in Fig. 15 
indicates twice the total experimental uncertainty of AQ^d^g ata , corresponding to plus and 
minus one standard deviation. 

It is seen that A<5j^ ad ^ 1 g ata can be reduced to zero for Aqcd ~ 1.3 GeV. As Aqcd is increased, 
more soft gluons are produced, increasing the probability for multiple color reconnections in an 
event. In events with multiple reconnections, color strings can reconnect the isolated gluonic 
string segment illustrated in Fig. lb back with the rest of the event, spoiling the rapidity gap. 
From Fig. 15 it is also seen that A<5^ 1 ^ 1 < g ata can be reduced to near-zero for large values of Q , 
e.g. Qo > 4 GeV/c 2 . As Qq is increased, fewer soft gluons are available, effectively decreasing 
the reconnection probability. In this sense, an increase in the value of Qo is analogous to a 
reduction in the value of the parameter R discussed above. We note that the values of Aqcd 
and Qo required to reduce AQ^^^ to zero or near-zero represent large excursions from their 
standard values. Fig. 15 suggests it is unlikely that AQj^ a ^g ata can be reduced to zero or 
near-zero through variation of a or a q . 

Setting Qo to 3.5 GeV/c 2 with the other parameters at their standard values, the Rathsman- 
CR model predicts a mean charged multiplicity in inclusive Z° events of (n ch .) =20.2, smaller 
than the experimental result of 21.15 ± 0.29 mentioned in Sect. 4. Mean multiplicity in the 
Lund hadronization model is primarily controlled by the parameters a and b. Therefore, having 
set Qo = 3.5 GeV/c 2 , we varied the b parameter 7 to reproduce the measured result for (n c h.). 
To increase the prediction for (n ch .), b needs to be decreased. As b decreases, AQ{^~ n g ata also 
tends to become smaller (see Fig. 15d). By iterating the adjustment of Q and b, it therefore 
proved possible to simultaneously obtain A(5^ 7 n ' g ata ~ and (n ch ) ^ 21. 15. The result we find 
for the two parameters is Q = 5.5 GeV/c 2 and b = 0.27 GeV~ 2 , corresponding to AQ^!^g ta = 
6.7 x 10~ 7 . We refer to the Rathsman-CR model with these adjusted parameters as the "re- 
tuned" Rathsman-CR model. 

A one standard deviation limit was evaluated for the re-tuned parameters by adjusting Qo 
and b so they yielded the correct result for (n^ ) and agreement with the one standard deviation 
upper limit for AQ|^7 n d g ata shown in Fig. 15: the result is Q = 3.7 GeV/c 2 and b = 0.35 GeV~ 2 . 
A two standard deviation limit was evaluated in an analogous manner, based on twice the total 
uncertainty of AQ^rJ ata : the result is Q = 3.2 GeV/c 2 and 6 = 0.38 GeV~ 2 . Finally, Q and 
b were adjusted to yield AQ|^j^ 1 g ata ~ and a value for (n c h.) equal to the LEP-averaged result 
plus its one standard deviation uncertainty, specifically (n c h.) = 21.44, see above and Sect. 4. 
The result is Q = 5.2 GeV/c 2 and b = 0.25 GeV" 2 . 

We examined the description of the re-tuned Rathsman-CR model for the inclusive Z° 
measurements in Figs. 2-5. The total x 2 f° r the 81 bins of data was found to be 1117.7, 
much larger than the result x 2 = 243.5 presented in Sect. 4 for the standard version of the 
Rathsman-CR model (see Table 3). Using the one and two standard deviation limits for the 
re-tuned parameters, given above, the corresponding x 2 are 435.1 and 327.2 respectively, still 
significantly larger than the x 2 of the standard version. The x 2 result for the parameters tuned 
to yield (n ch .) = 21.44 is 785.2. 

7 The a and b parameters are highly correlated with respect to the model predictions for (n c h.); therefore we 
consider variation of the b parameter alone, not both a and b. 
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Figure 15: Results for the difference between the Rathsman-CR Monte Carlo prediction and 
the experimental result for the Qi ea din g = bin in Fig. 13b, AQ lead T 1 g ata , as the principal pa- 
rameters of the model are changed with the other parameters at their default settings. The 
solid dots indicate the standard values of the parameters. The horizontal error ranges show 
the uncertainties attributed to the parameter values in [23]. For the parameters Aqcd and b, 
these uncertainties are too small to be visible. Note an uncertainty is not evaluated for the a 
parameter. The width of the shaded band centered on A(5^ ad ^ g ata = equals twice the total 
experimental uncertainty of AQ lead ^ g ata . 
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The x 2 results for the re-tuned Rathsman-CR model are listed in the bottom portion of 
Table 3. The deviations of the re-tuned model from the measured distributions are shown by 
the dotted curves in part (b) of Figs. 2-5. 

We attempted to follow an analogous procedure to that described above for Q to ad- 
just Aqcd- With Aqcd set to 1.3 GeV and the other parameters at their standard values, the 
mean charged multiplicity of inclusive Z° events is 26.4. To reduce this to 21.15, we increased 
the value of b. As b increases, AQj^^g ata becomes larger, however (Fig. 15d), and we could 
not find a solution which yielded both AQ^^f 1 ^ ~ and (n c h.) ~ 21.15. The closest solution 
we found, defined by the set of parameters which provided the correct inclusive multiplicity 
and a minimal result for AQ^J^, was A QCD = 1.3 GeV and 6 = 4.9 GeV -2 , which yielded 
A <3iid7ng ata = 0-012, with x 2 = 2-1 x 10 4 for the data of Figs. 2-5. 

Last, motivated by the observation that larger values of Aqcd and Q both reduce A<5^ ad ^ 1 ' g ata 
while having opposite effects on (n c h.), we increased both A QCD and Q , with the other param- 
eters at their standard values, to search for a solution with A(5^ n g ata «0 and (n ch .) ss 21.15. 
Specifically, we systematically increased the value of Aqcd and then performed a scan to de- 
termine the value of Qo which yielded the correct result for (n c h.). The parameter set with 
the minimal result for AQ|^7 n d g ata was A QCD = 0.6 GeV and Q = 4.9 GeV/c 2 , which yielded 
A <3iidr n d g ata = 0-0° 4 > with X 2 = 1606 for the data of Figs. 2-5. 

We conclude it is unlikely that the gluon jet results of Sect. 7.3 can be reproduced by the 
Rathsman-CR model through variation of Aq C d, similar to our observation above for a and a q . 

Thus, the only mechanism we found to adjust the parameters of the Rathsman-CR model 
to simultaneously describe our data on rapidity gaps in gluon jets and (n ch ) in inclusive Z° 
decays was to increase Qo to values in the range from about 3.3 to 5.5 GeV/c 2 , much larger 
than the values of 1-2 GeV/c 2 normally attributed to this parameter. These large values of Q 
resulted in a significant degradation of the model's description of inclusive Z° events, however, 
as discussed above. We conclude it is unlikely that this model can simultaneously provide a 
satisfactory description of the data in Sects. 4 and 7.3 using its standard value for the strength 
of color reconnection, Rq = 0.1. Thus, our results provide compelling evidence to disfavor color 
reconnection as it is currently implemented by this model. 



The Ariadne-CR model: The following parameters of the Ariadne-CR model were varied 
to determine their influence on AQj^Tg ata : 

• Aqcd, given by PARA(l); 

• Pr,min., the minimum transverse momentum of a gluon with respect to the dipole which 
emits it, given by PARA (3); 

• a, b and a q , given by PARJ(41), PARJ(42) and PARJ(21) as for the Rathsman-CR model. 

The PARA references are the names of the parameters in Ariadne, see Table 2. Analogous to our 
treatment of the Rathsman-CR model, we do not include the color suppression factor PARA(26) 
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(see Sect. 3) in this list. Varying PARA(26) to reproduce the result for the Q leading = bin in 
Fig. 13b, as well as the result for the Qi ea ding — bin plus its one standard deviation total 
uncertainty, yields PARA(26) =96 and 41, respectively, much larger than the standard value 
PARA(26) = 9. Note that large values of PARA(26) correspond to the limit of large Nc in 
which the probability for color reconnection becomes negligible. 

The results are shown in Fig. 16. The standard parameter values are indicated by solid 
dots. Their uncertainties as given in Table 2 are too small to be visible in the figure. Note 
an uncertainty was not evaluated for the a or a q parameters. Also note the Ariadne computer 
program requires Aqcd < PT,mm.- For this reason, the results for Aqcd are shown up to 0.70 GeV 
only, which is the standard value of pr,min.- 

The results of Fig. 16 are similar to those of Fig. 15, i.e. AQ^^n^ approaches zero as 
the parton shower cutoff pr.mm. is increased from its standard value, while it exhibits the same 
behavior shown in Fig. 15 as a, 6 and a q are varied. 

Setting p T ,mm. = 2 GeV/c so that AQ^ ad ^g ata ^ (see Fig. 16b), with the other parameters 
at their standard values, (n cn .) in inclusive Z° decays is predicted to be 20.0. Through iterative 
adjustment of p T , min . and 6, we found AQ^ c d r d g ata « (specifically, Ag^ a c d 7 n d g ata = -4.1 x lO -5 ) 
and (n ch .) « 21.15 for p T ,min. =4.7 GeV/c and 6 = 0.17 GeV -2 . We refer to the Ariadne-CR 
model with these adjusted parameters as the "re-tuned" Ariadne-CR model. One and two 
standard deviation limits were evaluated for the re-tuned parameters in the same manner as 
described above for the re-tuned Rathsman-CR model; the results are PT,min. = 2.0 GeV/c and 
6 = 0.35 GeV -2 , and PT,min. = 1-5 GeV/c and 6 = 0.42 GeV~ 2 , respectively. Similarly, we tuned 
Pr,min. and 6 to yield AQ^7 n g ata ~ and (n ch ) = 21.44, analogous to the procedure in Sect. 7.4: 
this yielded p TMn . = 3.5 GeV/c and 6 = 0.20 GeV -2 . 

Using the re-tuned parameters to determine the predictions of the model for the data in 
Figs. 2-5 resulted in a total x 2 of 3019.3 for those distributions, compared to the result x 2 = 32.4 
presented in Sect. 4 for the standard version of Ariadne-CR (Table 3). For the one and two 
standard deviation re-tuned parameters, the corresponding results are x 2 = 333.1 and 132.6, 
respectively. Thus the description provided by the re-tuned model is much worse than that 
provided by the standard version, even considering the uncertainties of the re-tuned parameter 
set. For the parameters tuned to yield (n ch .) = 21.44, the corresponding x 2 is 1254.6. 

The x 2 results for the re-tuned Ariadne-CR model are listed in the bottom portion of 
Table 3. The deviations of the re-tuned model from the measured distributions are shown by 
the dotted curves in part (c) of Figs. 2-5. 

Analogous to the procedure we followed for the Rathsman-CR model, we also attempted 
to simultaneously describe our gluon jet data and (n c h.) in inclusive Z° decays by varying 
both Aqcd and pT,mm. with the other parameters at their standard values. The solution which 
yielded the minimal result for A(5^ d ^ 1 g ata while correctly describing (n c h.) was Aq CD = 0.5 GeV 
and p T ,min. = 2.4 GeV/c, which yielded AQ%£j^** = 0.0041, with X 2 = 2095 for the data of 
Figs. 2-5. We did not attempt an adjustment of the two parameters Aq CD and 6 as we did for 
the Rathsman-CR model since the value of Aq CD is constrained by PT,mm. as explained above. 

Thus the only manner we found to adjust the parameters of the Ariadne-CR model to 
describe our data on gluon jets with a rapidity gap and at the same time yield (n cn .) ~ 21.15 



34 



0.04 

-2 

- = 0.02 

al 

C 

<l 

-0.02 



1 1 1 1 1 1 1 1 1 1 1 1 


1 1 


f(a) 




i< — limit from 




P T ,min 

1 




1 
1 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 







A QCD (GeV) 



0.04 

5 M 

-=0.02 

§ 41 

0> 

<l 

-0.02 



« 2f>, 
{j— 



6 o 

< 

-0.02 



1 1 1 1 1 1 1 1 1 1 1 1 

f(c) 


i i 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 


12 
a 


3 


i i i 1 i i i 1 i i i 1 i i i 1 
f(e) 


i i i 






i i i 1 i i i 1 i i i 1 i i i 1 


i i i 



0.2 0.4 0.6 0.8 
a q (GeV/c) 



0.04 
0.02 


-0.02 



0.04 
0.02 


-0.02 





i i i i i i i i i i 


r(b) 












i i i 1 i i i 1 i i i 1 





2 4 6 




P T , ml „ (^V/c) 


1 1 1 1 1 


i i i i | i i i i | i i i i | i i i i 


i i i i i i i 1 1 i i i i i i 1 1 i i i i i i i 







2 3 2 

b (GeV 2 ) 



Ariadne- CR 
versus data 



Figure 16: Results for the difference between the Ariadne-CR Monte Carlo prediction and the 
experimental result for the Qieading = bin in Fig. 13b, AQ lead r g ata J as the principal parameters 
of the model are changed with the other parameters at their default settings. The solid dots 
indicate the standard values of the parameters. The uncertainties attributed to the parameter 
values in Table 2 are too small to be visible. Note an uncertainty is not evaluated for the a 
or Oq parameters. The width of the shaded band centered on AQ^^g ata = equals twice the 
total experimental uncertainty of A<5{^^t ata - 
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was to increase the parton shower cutoff parameter pT,min. to values in the range from about 1.6 
to 4.7 GeV/c, significantly larger than the range of 0.6-0.7 GeV/c normally attributed to this 
parameter. This resulted in a significant degradation of the model's description of the global 
properties of inclusive Z° events as described above. Analogous to the Rathsman-CR model, 
we conclude it is unlikely that the Ariadne-CR model can simultaneously provide a satisfactory 
description of the data in both Sects. 4 and 7.3 using its standard value for the strength of 
color reconnection, and that our results provide compelling evidence to disfavor this model. 



8 Search for glueball-like resonances 

Besides providing a sensitive means to test models of color reconnection, gluon jets with a 
rapidity gap present an environment which may favor the production of glueballs, as discussed 
in [9]. If a hard, acollinear gluon in a e + e~ three-jet qqg event propagates a significant distance 
without radiating, a rapidity gap can form between the gluon jet and the rest of the event. 
This could enhance the probability for a color octet field to be created between the gluon and 
residual qq system, see Fig. 17a. This octet field is analogous to the field which is expected to 
connect two separating gluons produced in a color singlet state. This is in contrast to e.g. the 
Lund hadronization model, in which only color triplet fields are present, see Fig. 17b. Color 
octet fields provide a natural environment in which to create glueballs [40], through gg pair 
production from the vacuum, see Fig. 17c. 

QCD lattice calculations and other sources (see [41] and references therein) suggest that the 
mass of the lightest glueball state should lie in the general range from about 1 to 2 GeV/c 2 . 
One of the main candidates for the lightest glueball is the / (1500) [42], with a resonance width 
of 0.11 GeV/c 2 . The principal charged particle decay modes of the / (1500) are 7r + 7r~ and 
7T + 7r~7r + 7r~. 

Following the suggestions in [9] , we therefore examine invariant mass spectra in the leading 
part of our selected gluon jets. Since we are searching for anomalous resonant structure, the 
data are compared to the predictions of the models without color reconnection, i.e. Jetset, 
Ariadne and Herwig. These models do not contain glueballs. The distributions are examined 
at the detector level only and are normalized to the number of entries in the distributions. 
The bin widths are adjusted to reflect the mass resolution of the detector, estimated from the 
simulations. 

We begin by examining the total invariant mass of the leading part of the gluon jets, M leading . 
These masses are determined using both charged and neutral particles. Charged particles are 
assumed to be pions and neutral particles photons, as stated in Sect. 2. Since glueballs are 
electrically neutral, we select gluon jets with Qieading = 0, see Fig. 12a. Of the sample of 655 
events discussed in Sect. 6, this yields 250 gluon jets. The Mi ea di ng distribution of these jets is 
presented in Fig. 18a. The data are shown in comparison to the corresponding results of Jetset, 
Ariadne and Herwig. The models are seen to describe the data reasonably well. There is a 
general excess of the data above the Monte Carlo predictions for three bins in the mass range 
from 1.0 to 2.5 GeV/c 2 . The total x 2 values with respect to the data for these three bins are 
11.5 for Jetset, 7.3 for Herwig and 6.0 for Ariadne. Since this excess is only about two standard 
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(a) 



(b) 



/ggQg 
glueball 

(c) 



Figure 17: (a) Schematic illustration of a three-jet qqg event with a color octet field stretched 
between the gluon g and residual quark- ant iquark qq system. The double line with hatching, 
attached to the gluon, represents the octet field. The single lines connecting the octet field to 
the quark and antiquark represent color triplet fields, (b) Illustration of a qqg event in which 
color triplet fields connect the gluon directly with the quark and antiquark. (c) The octet field 
can be neutralized by the production of virtual gg color singlets from the vacuum, leading to 
the formation of glueballs. 



deviations of the statistical uncertainties above the predictions of Herwig and Ariadne 8 , it is 
not possible to obtain a definite conclusion concerning this discrepancy. 

Motivated by the charged particle decay modes of the /o(1500), mentioned above, we also 
examine the distributions of invariant mass of two oppositely charged particles in the leading 
part of the gluon jets, M^~ ding , and the corresponding distribution of four charged particles 
with total electric charge zero, M^~ d f n ~. The simulations predict that about 75% of the charged 
particles in the leading part of the selected gluon jets are pions. The distributions of M^^ ing 
and M^~ d f^ are presented in Figs. 18b and c. Since in this case the glueball candidate does 
not necessarily comprise the entire leading part of the gluon jet, there is no reason to constrain 
Qieading to be zero. Therefore, the M^~ ding and M^t" distributions are based on the entire 
sample of 655 gluon jets with a rapidity gap, not just the jets with (3i ea ding = 0. Again, the 
simulations are seen to describe the data reasonably well. The most significant excess of data 



above the Monte Carlo predictions occurs in the tail of the M } 



leading 



distribution, at mass 



values between 2.5 and 3.0 GeV/c 2 . This excess is about two standard deviations of the 
statistical uncertainties above the predictions of Ariadne and somewhat larger for the other 
models. Therefore, we do not observe clear evidence for anomalous resonant structure. Note 
the spike in the Herwig prediction for M^~ ding pa 0.77 GeV/c 2 in Fig. 18b is due to the p meson 
resonance which is too narrow in Herwig. Jetset exhibits a similar effect but at a less significant 



8 The Monte Carlo results are based on about twice as many events as the experimental distributions, see 
Sects. 2 and 3. 
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Figure 18: Distribution of jet invariant mass based on charged and neutral particles in the 
leading part of gluon jets. The distribution includes the effects of initial-state photon radiation 
and detector acceptance and resolution. The uncertainties are statistical only. The results are 
shown in comparison to the predictions of QCD Monte Carlo programs which include detector 
simulation and the same analysis procedures as are applied to the data. The hatched area 
shows the quark jet background evaluated using Herwig. 
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level. 

To better isolate a signal from a scalar particle such as the /o(1500), we also examined 
the cos 9* distribution of charged particle pairs in the leading part of the gluon jets, defined 
as follows. A "parent" momentum is defined by summing the momenta of two oppositely 
charged particles. 9* is the angle between the either of the two "decay" particles and the 
parent momentum, in the rest frame of the parent. The distribution of cos#* should be flat 
for a scalar particle but not necessarily for the combinatoric background. The measured cos 9* 
distribution was found to be well described by the Monte Carlo simulations. Therefore, we do 
not obtain any evidence for the anomalous production of scalar particles. 



9 Summary and conclusion 

A sample of 12 611 gluon jets with a mean energy of 22 GeV and estimated purity of 95% is 
identified in e + e~ hadronic Z° decay events using b quark jet tagging. The data were collected 
with the OPAL detector at LEP. A subsample of about 5% of these jets is selected which 
exhibit a rapidity gap, i.e. an absence of charged and neutral particles over a significant range 
of rapidity as illustrated in Fig. 7. After imposing the rapidity gap requirement, the estimated 
purity of the gluon jets is 86%. 

We examine the predictions of three models of color reconnection (CR): the Lonnblad 
model [25] (see also [28]) implemented in the Ariadne Monte Carlo, the Rathsman model [27] 
implemented in the Pythia Monte Carlo, and the color reconnection model in the Herwig Monte 
Carlo [18]. We refer to these as the Ariadne-CR, Rathsman-CR, and Herwig-CR models, re- 
spectively. Specifically, we examine the predictions of these models for the distributions of 
charged particle multiplicity and total electric charge in the leading part of the gluon jets, 
defined by charged and neutral particles beyond the gap. 

We find that the Rathsman-CR and Ariadne-CR models predict a large excess of gluon jets 
with a rapidity gap, for which the leading part of the jets is electrically neutral, compared to 
the corresponding models without color reconnection. In particular, these two models predict 
large spikes in the charged particle multiplicity distribution at even values of multiplicity. 
Thus our analysis is very sensitive to CR effects. We adjust the principal parameters of the 
two models to determine if they can be tuned to simultaneously provide a good description 
of our gluon jet measurements and the global properties of inclusive events in hadronic Z° 
decays. We find we can obtain a satisfactory description of the gluon jet data and the mean 
charged particle multiplicity (n c h.) in inclusive Z° events only for very large values of the parton 
shower cutoff parameters, Qo ~ 3.2-5.5 GeV/c 2 for the Rathsman-CR model or pr,min. ~ 1-5- 
4.7 GeV/c for the Ariadne-CR model, and that the overall description of global distributions 
in inclusive Z° events is then severely degraded. We conclude that it seems unlikely that either 
of these two models can be tuned to provide a satisfactory description of both our gluon jet 
data and the global properties of Z° events, using their standard values for the strength of 
color reconnection. We therefore conclude that color reconnection as currently implemented in 
these models is disfavored. Our conclusion for the Ariadne-CR model is consistent with our 
previous results [11]. Here, we present an even more sensitive study of color reconnection and 
systematically examine the effects of parameter variation on the model's predictions. 
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The Herwig-CR model also predicts a significant excess of events at even values of charged 
particle multiplicity in the leading part of gluon jets, compared to the corresponding model 
without color reconnection, cf. Fig. 14a. These excesses are much less prominent than for the 
Rathsman-CR and Ariadne-CR models, however, and are not clearly visible once the effects 
of finite detector resolution are incorporated, cf. Fig. 11. Therefore, we are unable to obtain 
a definite conclusion concerning this model. The data are nonetheless better described by the 
version of Herwig without color reconnection. 

Our study is also potentially sensitive to the presence of color singlet, electrically neutral 
objects such as glueballs, see [9]. We therefore examine the total invariant mass distribution of 
the leading part of gluon jets, using events in which the leading system is electrically neutral. 
We also examine the invariant mass distributions of two oppositely charged particles, and of four 
charged particles with total electric charge zero, in the leading part of our sample of gluon jets. 
We do not observe any evidence for anomalous features in the data, including the production 
of glueball-like objects. 
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